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ABSTRACT 

We present a new procedure to fit non-axisymmetric flow patterns to 2-D velocity maps of spiral 
galaxies. We concentrate on flows caused by bar-like or oval distortions to the total potential that may 
arise either from a non-axially symmetric halo or a bar in the luminous disk. We apply our method 
to high-quality CO and Ha data for the nearby, low-mass spiral NGC 2976 previously obtained by 
Simon et al., and find that a bar- like model fits the data at least as well as their model with large 
radial flows. We find supporting evidence for the existence of a bar in the baryonic disk. Our model 
suggests that the azimuthally averaged central attraction in the inner part of this galaxy is larger 
than estimated by these authors. It is likely that the disk is also more massive, which will limit the 
increase to the allowed dark halo density. Allowance for bar-like distortions in other galaxies may 
either increase or decrease the estimated central attraction. 

Subject headings: galaxies: kinematics and dynamics — galaxies: structure — galaxies: individual 
(NGC 2976) — galaxies: spiral — dark matter 



1. INTRODUCTION 

One of the first steps toward understanding the for- 
mation and evolution of galaxies is a determination of 
the radial distribution of mass within a representative 
sample of systems. The rotational balance of stars and 
gas in the plane of a disk galaxy offers a powerful probe 
of its mass distributi on, and has been widely exploited 
(|Sofue fc Rubinll200lb . 

When the motions of these tracers are consistent with 
small departures from circular orbits, the determination 
of the rotation curve (more precisely, the circular orbital 
speed profile) is s traightforwar d. However, it has long 
been known (e.g. iBosmal fl978T ) that large non-circular 
motions driven by bar-like or oval distortions, warps, 
or lopsidedness are common features in galaxy velocity 
maps, which complicate the determination of the radial 
mass profile. Yet the observed flow pattern contains a 
great deal of information about the mass distribution, 
which we wish to extract from the data. 

Since galaxies with closely flat, nearly axisymmetric 
disks are the exception, it is desirable to be able to esti- 
mate a mass profile in the more general cases. A number 
of techniques, which we review in ij2j already exist for this 
purpose. A procedure for dealing with a warped di sk has 
been successfully developed (e.g. lBegemanMl987D from 
the fi rst simple tilted ring analyses fe.g. iRogstad et ahl 
11974ft . and is now widely used. 

Non-axisymmetric distortions to the planar flow can 
always be described by an harmon i c anal ysis. But the ap- 
proach pioneered by iFranx et al.l (|1994l ) for interpreting 
the resulting coefficients embodies epicycle theory, which 
is valid only for small departures from circular orbits and 
may give misleading results if the observed non-circular 
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motions are not small compared with the circular orbital 
speed. A number of authors (see for references) ap- 
pear to find significant radial flows with speeds that rival 
the inferred mean orbital motion. Such flows violate the 
assumption of small departures from circular motion, are 
physically not well motivated, and the results are hard 
to interpret. 

We therefore propose here a new technique for fitting 
a general non-axisymmetric model to the velocity field 
of a galaxy that allows for large non-circular motions. 
We develop and apply the method specifically for the 
case of bar-like or oval distortions, but the procedure is 
readily generalized for potentials having other azimuthal 
periodicities. 

Our simple kinematic model, which we describe in $3l 
yields both the mean orbital speed and the amplitudes 
of the non-circular streaming velocities. It is successful 
because (1) we invoke a straight bar-like distortion to the 
potential, (2) we do not need to assume small departures 
from circular motion, and (3) we fit to the entire velocity 
field at once. 

We apply our method (£|4|) to the high-quality veloc- 
ity inap^_of_^GC 2976 that were previously presented 
bv lSimon et~aT1 (|2003l hereafter SBLB), and find that it 
suggests a significantly different radial mass profile from 
that deduced by those authors. We show (SJSJ) the reason 
for this difference, and argue that a bisymmetric distor- 
tion is both a more reasonable physical model, and that 
it is supported by independent evidence of a bar in this 
galaxy. 

2. MODELING NON-AXISYMMETRIC FLOWS 

2.1. Mathematical preliminaries 

The velocity of a star or an element of gas in the plane 
of the disk of a galaxy generally has two components 
at each point: tangential, Vt, and radial, V r , relative to 
any arbitrary center, most conveniently the kinematic 
center. Without approximation, each component can be 
expressed as a Fourier series around a circle of radius r 
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in the disk plane: 

oo 

V t (r, 9) = V t (r) + V rn.t{r) cos [mO + m>t (r)] (1) 

m—1 

and 

oo 

V r {r, 9) = V r (r) + £ V m>r (r) cos [mfl + 0m,r(r)] , (2) 

m— 1 

where the coefficients, V m ,t and V^r, and phases relative 
to some convenient axis, 6 m j and m ,r, are all functions 
of r. The quantity Vt(r) is the mean streaming speed of 
the stars or gas about the center; throughout, we refer to 
this quantity as the mean orbital speed. The axisymmet- 
ric term of the radial motion, V r (r), represents a mean 
inflow or outflow in the disk plane, which gives rise t o a 
"continuity problem" if it is large (| Simon et~aT1[2005h . 

Galaxies are observed in projection, with inclination i, 
about a major axis, which we choose to define as 9 = 
in the above expansions. The line-of-sight velocity is 
the sum of the projected azimuthal and radial velocities: 
V hs — V sys + sin i(Vt cos 9 + V r sm9), where V^ ys is the 
systemic velocity of the galaxy. In terms of our Fourier 
series, 

^obs — ^sys 

{oo 
Vt COS 9 + Vrn.t cos 9 cos [mO + 9 mt t] 
m=l 

oo ^ 

+ V r sin 9 + Vm.r sin 9 cos [m9 + 9 rn , r \ > . (3) 

m—l ) 

Using standard trigonometric relations, this expression 
can be rewritten as 

^obs ^sys T ~> n 

: — = V t cos 9 

smi 

+ £^{cos[(m + l)0 + m , t ] 

m—l 

+ cos [(m- 1)6* + 9 mjt ]}+ V r sin9 
+ £ ^{sin[(m + l)0 + m , r ] 

m—l 

- sin[(m- l)9 + 9 m . r }}. (4) 

As is well known (e.g.[C anzian 1993; Schoenmakers et all 
[19971 iCanzian fc Allenlll997t iFridman et al.1 12001ft . pro- 
jection therefore causes velocity distortions with intrinsic 
sectoral harmonic m to give rise to azimuthal variations 
of orders m! = m ± 1 in the corresponding line-of-sight 
velocities. Thus intrinsic distortions at two different sec- 
toral harmonics give rise to projected velocity features 
of the same angular periodicity in the data, complicat- 
ing the determination of all coefficients in the expansion. 

2.2. Previous approaches 

The principal scientific objective of most spectroscopic 
observations of disk galaxies is to extract the func- 
tion Vj(r), which should be a good approximation to 
the circular orbital speed if all other coefficients on 
the right-hand side of eq. [3] are small. With a sin- 
gle slit spectrum along the major axis of the galaxy, 



one generally sets V b s = V aya + Vt(r) sin i, implic- 
itly assuming all other terms to be negligible. In this 
case, the inclination must be determined from other 
data (e.g. photometry). Differences larger than mea- 
surement errors between the approaching and reced- 
ing sides flag the existence of non-circular motions, 
but measurements along a single axis do not yield 
enough information to determine any other coefficient. 
This and other uncert ainties inherent in such deduc- 
tions are well- r ehear s ed (Ivan den Bosch fe Swatersll2001l; 
de Blok et al.l l2003t ISwaters et al.l l2003al iRhee et all 
20041 : iSne kkcns et al '. 120051 : lHavashi fc Navarroll2007h . A 
two-dimensional velocity map, on the other hand, pro- 
vides much more information. 

Software packages, such as rotcur (Bcgcman 1983), al- 
low one to fit the velocity field with a single velocity 
function Vt(r) in a set of annuli whose centers, posi- 
tion angles (PAs) and inclinations are allowed, if desired, 
to vary with radius. This package is ideal for the pur- 
pose for which it was designed: to determine the mean 
orbital speed even when the plane of the disk may be 
warped. It works well when non-circular motions are 
small, but yields spurious variations of the parameters 
when the underlying flow contains non-axisymmetric, es- 
p ecially bisymmetric, distort ions. 

iBarnes fc Sellwoodl (|2003h adopted a different ap- 
proach. They assumed the plane of the disk to be flat, 
and determined the rotation center, inclination, and PA 
by fitting a non-parametric circular flow pattern to the 
entire velocity map. Their method averages over veloc- 
ity distortions caused by spiral arms, for example, but 
again may yield spurious projection angles and mean or- 
bital speeds if there is a bar-like or oval distortion to 
the velocity field over a wide radial range. The rotcur- 
shape program, recently added to the NEMO (|Teubenl 
1995) package, suffers from the same drawback because 
it also assumes a fiat, axisymmetric disk. Furthermore, 
it fits multiple parametric componen ts to a velocity field 
and th us has less flexibility than the IBarnes fc Sellwoo"dl 
<|2003T) techniq ue. 

iFranx et all (|1994h and ISchoenmakers et all (|1997l ) 
pioneered efforts to measure and interpret the non- 
axisymmetric coefficients that describe an observed ve- 
locity field, and expansions up to order m' ~ 3 
are now routi nely c a rried out (e . g. IWong et alj [20041 
Chemin et ail l2006t ISimon et all 120051: iGentile et all 
Their approach assumes departures from circu- 
lar motion to be small so that the radial and tangential 
perturbations for any sectoral harmonic can be related 
through epicycle theory. The technique is therefore ap- 
propriate only when all fitted coefficients are small and 
the mean orbital speed is close to the circular orbital 
speed th at balances the azim uthally averaged central at- 
traction. I Wong et al.1 (|2004l ) present an extensive discus- 
sion of this technique and conclude that it is difficult to 
work backwards from the derived Fourier coefficients to 
d istinguish between diff e rent ph ysic al models. 

ISwaters et al.l (|2003bft . lSBLBl and IGentile et al.l (|2007l ) 
report velocity fields for nearby galaxies that show non- 
circular motions whose amplitude rivals the mean orbital 
sjpeed at small r. Swaters et al. (2003b) note that their 
Vt (r) model fails to reproduce the inner disk kinematics 
of their target. They correct Vj(r) for an isotropic veloc- 
ity dispersion of 8 km/s, but do not attempt to model 
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the isovelocity twists in their Ha velocity field. ISBLBl fit 
the simplest acceptable model to their data: an_axisym- 
mctric flow with just two non-zero coefficients Vt(r) and 
V r (r). They favor this model over a bar- like distortion 
partly because the galaxy is not obviously barred, and 
partly because they find that the rn! = 3 components 
are scarcely larger than the noise (see also Q. The ad- 
dition of the radial velocity term V r allows a more com- 
plicated flow pattern to be fitted with an axisymmet- 
ric mo del, which significan tly improves the fit to their 
data. iGentile et alj (|2007T ) do detect a radial m! = 3 
component in addition to a strong radial m' = 1 term 
in the kinematics that they report, which they conclude 
"are consistent with an inner bar of several hundreds of 
pc and accretion of material in the outer regions" . De- 
spite finding large non-circular motions, the authors of 
all three studies nonetheless adopted their derived mean 
orbital speed as the "rotation curve" of the galaxy, which 
they assume results from centrifugal balance with the az- 
imuthally averaged mass distribution. 

These deductions are suspect, however. As we show 
below (3S]l), a bisymmetric distortion to the flow pattern 
may not give rise to a large m' — 3 term in the velocity 
field, and the smallness of these terms does not establish 
the absence of a strong bisymmetric distortion. Further, 
associating the V t term with the rotation curve is valid 
only if all departures from circular motion are small, yet 
they had found non-circular velocity components almost 
as large as the mean orbital speed over a significant radial 
range. 

Early work on modeling gas flows in barred galaxies 
is reviewed i n Sellwood fc Wilki n son (1993; s e e thei r 
sect ion 6.7). iWeiner etlll (|200lD . iKranz et all (|2003f) . 
and lPerez et alj ()2004f ) attempt to build a self-consistent 
fluid-dynamical model of the non-axisymmetric flow pat- 
tern. They estimate the non-axisymmetric part of the 
mass distribution from photometry and try to match the 
observed flow to hydrodynamic simulations to determine 
the amplitude of the non-axisymmetric components of 
the potential. The objective of this, altogether more am- 
bitious, approach is to determine the separate contribu- 
tions of the luminous and dark matter to the potential. 
Here, our objective is more modest: to estimate the mean 
orbital speed from a velocity map that may possibly be 
strongly non-axisymmetric. Thus their attempt to sepa- 
rate the baryonic from dark matter contributions seems 
needlessly laborious for our more limited purpose. 
3. A NEW APPROACH 
3.1. A bar-like distortion 

Our objective is to model non-circular motions in a 2- 
D velocity map. Since we do not wish to assume that 
non-circular motions are small, we refrain from adopting 
the epicycle approximation. However, we do make the 
following assumptions: 

• The non-circular motions in the flow stem from a 
bar-like or oval distortion to an axisymmetric po- 
tential. We suppose these motions to be caused 
by either a non-axially symmetric halo in the dark 
matter or by a bar in the mass distribution of the 
baryons. 

• A strong bisymmetric distortion to the poten- 
tial, even one that is exactly described by a 




Fig. 1. — Parameter definitions and flow pattern in the disk plane 
for the bisymmetric model (eq. [5}. The solid circle represents the 
largest r included in the model, and the horizontal dash-dotted line 
is the major axis of the disk defined by the sky plane. The long- 
dashed line is the major axis of the bisymmetric distortion, at angle 
4>h from the major axis. The diamond denotes the location (x e , y e ) 
of a datapoint D n , a distance r n from the kinematic center (x c , y c ) 
and at PAs 9^ from the bisymmetric distortion axis and 9 from the 
major axis. The dotted circle shows the circular orbit of radius 
r n in the disk, and the axisymmetric model component Vt(r n ) is 
in the counter-clockwise direction. The extrcma of components 
V'2,t(i'n) and V2,r(rn) are indicated by solid and dashed arrows, 
respectively, and large dots at the same distance from (x c ,y c ) as 
each set of arrows denote PAs where the corresponding component 
passes through zero. 

cos(26*) angular dependence, can give rise to more 
complicated motions of the stars and gas (e.g. 
ISellwood fc Wilkinson] Il993f) . In particular, the 
flow may contain higher even harmonics. Never- 
theless, the m = 2 terms will always be the largest, 
and we therefore begin by neglecting higher har- 
monics. 

• We assume the bar-like distortion drives non- 
circular motions about a fixed axis in the disk 
plane. In a steady bar-like flow, the perturbed 
parts of the azimuthal and radial velocities must 
be exactly out of phase with each other. That 
is, the azimuthal streaming speed is smallest on 
the bar major axis and greatest on its minor axis, 
while radial motions are zero in these directions 
and pe ak, with alternating signs, a t intermediate 
angles (|Sellwood fc Wilkinsonill993f) . 

• We must assume the disk to be flat, because we 
require the predicted V b s from eq. [3] to have the 
same inclination at all r. This assumption is ap- 
propriate for spiral galaxy velocity fields measured 
within the op tical radius, where warps are rare (e.g. 
iBriggd fT990f ) . and is therefore well-suited to inter- 
preting kinematics derived from Ha, CO or stellar 
spectroscopy. The technique presented here should 
therefore not be applied to the outer parts of H I 
velocity fields, w hich typically extend w ell into the 
warp region (e.g. iBroeils fc Rhelfl99l . 
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A model based on these assumptions predicts the ob- 
served velocity at some general point in the map to be 
given by eq. [3J with the m — 2 terms as the only non- 
axisymmetric terms: 



Knodel = Kys + sin i [ Vt cos V 2 , t cos(26> & ) cos 9 
— V2, r sin(2^b) sin 9 1 



(5) 

The geometry in the disk plane is sketched in Fig. [1] As 
above, 9 is the angle in the disk plane relative to the 
projected major axis, which is marked by the horizon- 
tal dash-dotted line. The major axis of the bisymmetric 
distortion (or bar for short) lies at angle to the pro- 
jected major axis; thus angles relative to the bar axis are 
9b = 9 — 4>b- We have chosen the phases of the V^t and 
V2. r terms such that both are negative at 9b = and 
they vary with angle to the bar as the cosine and sine 
of 29b respectively. Comparing eqs. [HJandEJ we see that 
#2,t = 7r — 2(f>b and $2, r = 7r /2 — 2<\>b- 

The amplitudes of the tangential and radial compo- 
nents of the non-circular flow, V2,t(f) and V2. r (r) respec- 
tively, are both functions of radius. While it is possible 
to relate the separate amplitudes for a given potential, 
we do not attempt to model the mass distribution that 
creates the flow and therefore allow them both to vary 
independently. 

Circles in the disk plane project to ellipses on the sky, 
with ellipticity td given by 1 — ea = cosi, with a com- 
mon kinematic center, (x c ,y c ). We use primes to denote 
projected angles onto the sky plane. Thus the projected 
PA 3 of the disk major axis is <j>' d , while 4>' b is the PA of 
the bar major axis in the sky plane. These angles are 
related by 



4>'b = 4>'d + arctan(tan0f,cosi) 



(6) 



In addition to the three velocity functions Vt(r), p2,t(r) 
and V^rfr)) the model is therefore described by the pa- 
rameters (x c , y c , V sys , ed, 4»' d , 4>' b )- We refer to the model 
described by eq. [5] as the bisymmetric model. 

3.2. Other possible models 

Other models for the flow pattern could readily be de- 
rived from eq. [3j 

In particular, and solely to facilitate comparison with 
other work, we also fit a purely axisymmetric model with 
the coefficents of all m > terms set to zero, but retain 
the V r term. There is no undetermined phase angle for 
this intrinsically axisymmetric model and the predicted 
velocity is simply 



model 



= V^y S -h sin z [V t cos + V r sin 6 ] 



(7) 



The coefficient V r corresponds to pure radial inflow or 
outflow. 4 We will refer to this as the radial model. 

Other, more complicated, models could also be fitted to 
data by retaining more terms as required, provided that 
an assumption is made about the radial dependence of 
the phases of the non- axisymmetric perturbations. The 
extension of these formulae to include other velocity field 
harmonics is straightforward, and we have tried doing so 
in some of our analyses (see $4]). 

3 All PAs are measured North — » East. 

4 It is not possible to distinguish between inflow and outflow in 
this model unless the side of the disk along the minor axis that is 
nearest to the observer can be determined independently. 



3.3. Discussion 

If the non-circular motions measured in some spirals 
do stem from bar-like or oval distortions, then the bisym- 
metric model has several advantages over both the radial 
model and also over epicyclic approaches for characteriz- 
ing these asymmetries. Since m' = 1 velocity field com- 
ponents can arise from either radial flows or a bisym- 
metric perturbation to the potential (eq. [4]), both the 
bisymmetric and radial models could produce tolerable 
fits to the same data. However, the bisymmetric model 
offers a more direct, unambiguous approach for identify- 
ing m = 2 distortions than does the radial model. 

Moreover, interpretations of velo city field harmon- 
ics that rely on epicycle theory (Franx et al.l 119941 : 
ISchoenmakers et al.lll997l : rCanzian fc Allenlll997t ) are ap- 
plicable only in the limit of a weak perturbation to the 
potential, whereas the components of our bisymmetric 
model are not restricted to mild distortions. We also 
note that since the bisymmetric model imposes a fixed 
4>' b on the non-circular flow pattern, it is not sensitive 
to m = 2 perturbations to the potential that are not in 
phase (such as spiral patterns). 

Finally, the bisymmetric technique is much simpler 
than fluid-dynamical modeling of the velocity field (see 
§2.2p , since it does not require (or yield) a model for the 
mass distribution. 

3.4. Fitting technique 

We attempt to fit the above kinematic models to obser- 
vational data by a n extension of the minimiz ation pro- 
cedure devised by iBarnes fc Sellwoo"dl {2003). In gen- 
eral, we need to determine the systemic velocity V sys , 
kinematic center (x c ,y c ), ellipticity e^, and disk PA cj> d , 
as well as M unknown radial functions V m ,t and V m ^ r 
(rn = and m > if desired) and the (fixed) PA(s), 9 m , 
of any non-axisymmetric distortions to the flow. 

We tabulate each of the M independent velocity pro- 
files at a set of concentric circular rings in the disk plane 
that project to ellipses on the sky with a common center, 
{ x c He)- Once these tabulated values are determined, we 
can construct a predicted Vmodci at any general point by 
interpolation. We difference our model from the data, 
which consist of N line-of sight velocity measurements 
{D n } with uncertainties {cr n }, and adjust the model pa- 
rameters to determine the minimum m in °^ ^ ne stan- 
dard goodness-of-fit function Xt w hh v degrees of free- 
dom: 



(8) 



Here, the K elements of {T4} are the values of the tab- 
ulated velocity profiles in the model and the weights, 
Wk, n describe the interpolation from the tabulated Vk to 
Knodei (eq- El or [7|) at the position of the observed value 
D n . 

When Xr — Ar min' tne partial gradient of with re- 
spect to each Vj, where j labels each of the in turn, 
must satisfy 




dV 




Z)fe=l w k,nVk 



Wj 



0. (9) 
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Rearranging, we find 

k / n \ N 

E jata-U-ESffl., no) 

k=l \n=l / n=l 

resulting in a linear system of K equations for the K 
unknowns, {Vfc}- 

For a given set of attributes (x c , y c , V sys , £d, <f)' d , 8 m ) 
of the projected disk, we compute {14} by solving the 
linear system of eq.QIIl and use the resulting {14} values 
in eq.[5]to evaluate Xi- The best fitting model is found by 
minimizing eq. [5] over the parameters mentioned above, 
which necessitates recomputing {14} via eq. [10] at each 
iteration. Any convenient method may be used to search 
for the minimum' we use Powell's direction set method 
(iPress et al.lll992T). 

iBarnes fc Seilwool £2003) use this minimization strat- 
egy to extract only the mean orbi tal speed Vt(r) from 
Ha ve locity fields of spirals in the iPalunas fc Williams! 
(|2000l) sample. In our more general case, the M > 1 
model profiles are defined by distinct sets of K' M rings 
in the disk plane, and {14} contains all of the veloci- 
ties from these profiles: K = J^m^-m- 111 other words, 
adding a velocity profile (defined in K' rings) to a model 
increases the rank of the matrix in eq. [10] by K' . The 
radial model has M = 2, while M = 3 for the bisymmet- 
ric model. Further discussion of {14} and derivations of 
{wk,n} are given in the Appendix. 

4. VELOCITY FIELD MODELS OF NGC 2976 

To illustrate the technique, we fit our bisymmet- 
ric model to the obser ved h igh-quality velocity field of 
NGC 2976 reported by ISBLBL NGC 2976 is a nearby, 
low-mass Sc galaxy with i ~ 60°. We adopt a distance 
D = 3.5 6 Mpc, estimated from th e tip of the red giant 
branch (jKarachentsev et al.ll2002T ). and convert angular 
s cales to linear scales using 1" = 17.3 pc. 

ISBLBI present Ha and CO velocity fields of NGC 2976, 
with a spatial resolution 5 of ~ 5" (86 pc) and spectral 
resolutions of 13kms _1 and 2kms~ 1 , respectively. They 
find that the velocity field is not well-modeled by disk 
rotation alone. 

They report a detailed analysis of these kinematic data 
in which the projection geometry of their model rings is 
determined from optical and near-IR photometry. They 
conclude that there is no strong evidence for a bisym- 
metric distortion in this galaxy, since all m! > 1 compo- 
nents of the velocity field are consistent with noise. They 
find that a combination of rotation and pure radial flows 
provides an adequate fit. The amplitude of the inferred 
radial velocity profile rivals that of the rotational com- 
ponent for r < 500 pc: NGC 2976 thus exhibits some 
of the largest non-circular motions ever detected in a 
low-mass, rotationally-supported system. In their later 
paper, Simon et al. (2005) noted that finding large val- 
ues of the 14 term is a strong indication that a model 
with an axisymmetric radial flow is incorrect. They sug- 
gest that the non-circular motions in NGC 2976 stem 
from a t riaxial halo, but their use of epicycle theory re- 
lations _|SchoinmakersIi£]d] [1997]) is inappropriate be- 
cause V r (r) is not always small (see their fig. 9). We 

5 Throughout, we recompute the linear scales presented bv lSBLBl 
for consistency with our choice of D. 



also suspect that a bisymmetric distortion is responsible 
for the observed departures from a circular flow pattern. 

The Ha and CO velocity fields of NGC 2976 presented 
m fig. 4 of lSBLl were kindly made available to us by 
J. D. Simon. Following these authors, we analyse the 
kinematics of the two tracers together, since the data 
agree within their uncertainties. 

We fit the combined velocity field with our bisymmetric 
model to examine whethe r the departures from circular 
motion detected by SBLB stem from an m = 2 distortion 
to the potential. In order to demonstrate that our new 
technique (fJ3]) yi elds a similar kinematic model to the 
one obtained by ISBLBI . we apply our radial model to 
the same dataset, and compare values of the parameters 
we obtain as a consistency check on our method. For 
completeness, we also attempt to fit the data with a suite 
of other models including m = 0, m = 1 and m = 2 
distortions (see below). 

The observations presented by lSBLB j sample the veloc- 
ity field of NGC 2976 out tor - 130" (2.2 kpc) from its 
photometric center. We evaluate the velocity profiles in a 
maximum of K 1 = 26 rings, separated by 4" for r < 95" 
and by up to 10" farther out. Neither the bisymmetric 
model nor the radial model yielded reliable constraints 
on the non-circular components of the velocity field for 
r > 100". We therefore conclude that the outer part of 
NGC 2976 is adequately described by a simple circular 
flow, and fix the amplitudes of all coefficients but V t to 
zero beyond that radius. This reduces the rank of the 
matrix (eq. 1 1 0[) by a few. 

To check the validity of our planar disk assumption at 
the largest radii probed by the measurements, we com- 
pare the disk parameters derived from fits including dif- 
ferent numbers of outer rings. Specifically, each mini- 
mization uses the same ring radii, except that the outer- 
most ring included is varied in the range 80" < r max < 
135" and velocity measurements at radii beyond r max are 
ignored in the fit. Models with r max < 112" return iden- 
tical disk parameters within the uncertainties, but the 
optimal values of x c , y c and V sys change substantially 
when rings at larger r are added. We therefore restrict 
our fits to include only D n with r n < 112" in our final 
models, as the disk may be warped farther out. 6 

We make an allowance for ISM turbulence by redefin- 
ing {ct„} to be the sum in quadrature of the uncertain- 
ties in the emission line centroids and a contribution 
Aism = Skms^ 1 . We find that choosing values of Atsm 
in the range 3kms _1 < Atsm ^ 7kms _1 and varying 
the ring locations and sizes by 2-4" have little impact on 
our results. 

In addition to the bisymmetric and radial models of 
NGC 2976, we also fitted models including a lopsided 
(rn = 1) distortion. The optimal m = 1 model (includ- 
ing velocity profiles V t (r), V\ tt {r) and Vi >r (r); see eq. [3]) 
produced a much less satisfactory fit to the data than 
either the bisymmetric or the radial model. Adding a 
radial flow term V r (r) yielded optimal parameters iden- 
tical to those of the radial model, with the m = 1 com- 
ponents consistent with zero. The insignificance of a lop- 
sided component, which we conclude from these fits, is 

6 The disk geometry and kinematics of NGC 2976 at r > 1.5 kpc 
will be explored in detail using extant aperture synthesis H I maps 
of the system. 
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Fig. 2. — Kinematic models of the NGC 2976 velocity field. The panels in the top row show (a) the observed velocity field {D n } from 
ISBLBI . (6) the optimal radial and (c) the optimal bisymmetric models. The velocity fields are plotted on the same colorscale, shown in 

kms -1 to the right of the top row. The panels in the bottom row show residual maps {AV n } = {D n — EfcLi w k,nVk} for (d) the radial and 
(e) the bisymmetric models, with the colorscale in kms -1 for both shown to the right of that row. The model velocity fields and residuals 
have been rotated by — (ip' d + tt/2) about (x c ,y c ) from Table[T] The data in[2}i are rotated by the photometric value —(—37° + 7r/2) 
(intermediate to the two model values) about the photometric center 09 h 47 m 15! 3, 67° 55' 00'.' 4 ijSBLBT ). The orientation of {D„} (roughly 
correct for the models as well) and the map scale are at the bottom left. The unusual sha pe of c ontoured regions in the maps reflects the 
locations of the individual pointings used to construct the Ho velocity field of NGC 2976 ( S BLBI) , 



consistent with our result below that the kinematic and 
photometric centers of NG C 2976 are coincident within 
the errors (see also ISBL B). We also attempted to fit 
m = and m — 2 distortions to the data simultaneously. 
Since both cause m' = 1 periodicities in the linc-of-sight 
velocities, however, the resulting model had too much 
freedom and produced unphysically large variations in 
all the velocity profiles. 

4.1. Uncertainties 

The curvature of the Xr surface at the minimum im- 
plies small formal statistical errors on the best fitting 
model parameters because of the large numbers of data 
values. The x 2 m i n + 1-0/V contour on the Xt surface 
corresponds to variations SV < lkms -1 on the velocity 
profile points, which we regard as unrealistically small. 
We therefore use a bootstrap technique to derive more 
reasonable estimates of the scatter in the model param- 
eters about their optimal values. 

For each model, we generate a bootstrap sample of 
the data by adding randomly drawn residuals AV n = 

D n - EfcLi Wk,nVh from the distribution at x? = X?, m in 
to the optimal model velocity field. Since {AV n } is 
correlated over a characteristic scale corresponding to J 
adjacent datapoints, fully random selections do not re- 
produce the quasi-coherent residuals we observe in the 
data. We therefore select P = N / J values of AV n and 
add them to the model at P random locations drawn 
from {r n ,6?}; residuals at the remaining (1 — 1/J)N in 
{?"n,0} are fixed to the value of the nearest randomly 
drawn AV n . We find that J = 4 produces bootstrap 
residual maps with features on scales similar to those in 
{AV n } for the models of NGC 2976 in Fig. [2] (see below), 
but that there is little change in the derived uncertainties 
for 2 < J < 5. 

For both the bisymmetric and radial models, we there- 
fore construct bootstrap samples of the observed velocity 
field using J = 4. We repeat the minimization for each 
sample, substituting the bootstrap velocities for {D n } 
in eqs. |8] and HU We carry out this procedure 1000 
times, and adopt the standard deviation of each param- 
eter about its mean value from all the realizations as its 
1(7 uncertainty in the model of the measured velocities 



{D n }. 

4.2. Results 

Our final models of the SBLB Ha and CO velocity 
fields for NGC 2976 are shown in Fig. [5] The minimiza- 
tion results are given in Table [TJ and the corresponding 
velocity profiles are shown in Fig. (3[ 

The observed velocity field from lSBLBl is reproduced 
in Fig. [2]3, the best fitting radial and bisymmetric models 
are in Figs. [2^ and [2]c, and the residuals {Al^} are in 
Figs. EH and EJs. Both models reproduce the gross fea- 
tures of the observed velocity field, although the bisym- 
metric model exhibits a somewhat larger isovelocity con- 
tour "twist" along the kinematic minor axis (oriented 
vertically in Fig. [2]) than the radial model. The residual 
patterns in Figs. [2)5? andJSJs are very similar: AV n is cor- 
related on scales of 15 — 20" (250 — 350 pc) in the maps, 
which may reflect large-scale turbulence. The mean val- 
ues (|AV^|) in Table [T]are slightly lower for the bisym- 
metric model than for the radial one, as is also suggested 
by the colors in Figs. [2jf and[2]e. 

The values of (x c ,y c ) and V^ ys in the two models (Ta- 
ble [TJ) are identical within their uncertainties, while the 
radial model favors a larger and <p' d than the bisymmet- 
ric model at the 2<r level. Both sets of kinematic param- 
eters (x c , y c , €d, <t >ii) are consistent with the photometric 
values derived by ISBLBl . corroborating their conclusion 
that there is little evidence for an offset between them 
(see also The values of X?,min indicate that both 

models adequately describe {D n }, with x 2 ~ 1 per de- 
gree of freedom for the adopted Aism- 7 Even though the 
bisymmetric model has fewer v than the radial model, 
the difference between their x 2 m ; n is formally significant 
at the 12cr level. As with the unrealistically small model 
uncertainties implied by the X 2 m in + 1-0/z" contour on 
the Xt surface, however, a literal interpretation this dif- 
ference in goodness-of-fit is unwise. We thus conclude 
conservatively that the smaller x 2 m in an( i l° wer (|AV^|) 
of the bisymmetric over the radial model imply only a 

7 The optimal model parameters remain unchanged with choices 
of AigM in the range 3 kms - 1 < Atsm ^ 7kms _1 , b 1 ^ 2 
the corresponding fits varies from 2.7 < Xt min ~ 
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Fig. 4. — Difference V2,t{r) — V2,r(x) in the optimal bisymmetric 
model. The components are plotted separately in Fig. [3f>. 



Fig. 3. — Fitted velocity components for NGC 2976. (a) The components of the optimal radial model: Vt(r) is shown by the red circles 
and V r (r) is shown by the blue squares (eq. 0. (b) Velocity components from the optimal bisymmetric model: here Vj(r) is shown by the 
red circles, V2 t(r) by the green squares, and V2 ;r (r) by the blue triangles (eq.[5]l. 

with their results (the large Vp(r) fo r r < 10" in our 
radial model was also found by ISBLBl in their ringfit ve- 
locity_field decompositions; see their fig. 7a). We find 
that V r (r) is ~ 7kms _1 smaller at r < 30" t han the ra- 
dial velocity amplitudes presented by ISBLBl . This ~ 2a 
discrepancy results from our inclusion of a Aism term in 
{cr„} (eas. I8l— IlO]) . and disappears if we set Aism = in 
our model. Thus our analysis confirms the non-circular 
motions in NGC 2976 found by these authors. 

The bisymmetric model favors a strongly non- 
axisymmetric flow about an axis inclined 17° to the pro- 
jected major axis in the disk plane. The radial varia- 
tions and uncertainties of all three fitted velocity compo- 
nents are shown in Fig. [3j). The estimated uncertainties 
on the {Vfc} are larger than those in the radial model 
(Fig. (3]a), consistent with the larger scatter in the pro- 
file values from ring to ring. This is likely due to the 
extra velocity profile relative to the radial model, which 
gives the bisymmetric model increased flexibility to fit 
small-scale features. As in the radial model, we find 
significant non-circular motions, this time in the form 
of a bisymmetric flow pattern in the disk plane. The 
overall shape of the non-circular contributions V2,t(r) 
and V2,r(r) resembles that of V r (r) in the radial model: 
this is reasonable because both models must fit the 
m' = 1 variations of the velocity field. The difference 
V2,*(r) — V2 ;J .(r) between the bisymmetric components 
in Fig. [3] is plotted in Fig. [U There is marginal evi- 
dence that V 2 ,t(r) > V 2 , r (r) for 45 " < < 65", but 
elsewhere the two components have very similar ampli- 
tudes. Linear theory applied to a weak, stationary bar- 
like distortion produces V2,t (r) = V2, r {r) for a solid-body 
rotation velocity profile and t(r) = 1.5V2, r (0 f° r a 
flat one (|Sellwood k. Wilkinson! Il993l ). Although linear 
theory cannot be trusted for strong perturbations, it is 
somewhat reassuring that it predicts similar V 2 , t(r) and 
V2, r {r) for a rising V t (r). 

The most significant difference between the optimal 
radial and bisymmetric models is in the shape of Vt(r). 



(kpc) 




2.5 



Fig. 5. — Difference AVt(r) between the optimal Vt(r) from the 
bisymmetric model and the optimal Vt(r) from the radial model. 
The components are plotted separately in Fig. [3] 

marginally superior statistical fit to the data. 

The best fitting velocity field components from our ra- 
dial model are shown in Fig. [3)3. Despite significant dif- 
ferences between our minimization technique and that of 
SBLB, our measurements of Vt (r) and V r (r) agree well 
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Beyond the region affected by non-circular motions, r > 
80", Vt(r) is identical in the two models, as it must be, 
but large differences arise where non-circular motions are 
large. Fig. [5] shows the difference between Vt(r) from the 
bisymmetric model and that from the radial model: the 
former profile rises more steeply than the latter, and its 
amplitude is larger by ~ 15kms _1 for 15" < r < 50". 
We discuss the reason for these differences in the next 
section. 

5. DISCUSSION 
5.1. Mean streaming speed 

The large differences in the fitted Vt(r) for the bisym- 
metric and radial models over the inner part of NGC 2976 
(Fig. [5]) demand explanation. Fig. Oshows sky-plane pro- 
jections of the angular variations of the separate fitted 
velocity components at r = 20"_for both models. 

Fig. shows the projected V r (r) of the radial model 
(dash-dotted line), which shifts the peak in the projected 
model velocity (solid line) away from the kinematic ma- 
jor axis (9 = 0) and reproduces the iso- velocity "twist" in 
the observed velocity field (Fig. [2]). Since the projected 
V r must be zero at 9 = in this model (see eq. [7]), the 
projected Vt(r) (dashed line) must equal Knodei along the 
kinematic major axis. 

Fig. O) shows the corresponding case for the bisym- 
metric model where the non-axisymmetric terms allow 
Knodei to differ from Vt along the major axis (eq. [5]). 
The abscissae are marked as 9 along the bottom and as 
9t, along the top, which differ only slightly because the 
best fitting major axis of the bisymmetric distortion in 
NGC 2976 projects to a PA similar to the kinematic ma- 
jor axis (i.e. <f)' b ~ <fi' d ; Table [T]). The significant negative 
contribution of the projected V2,t(r) (dash-dotted line) 
to Vmodci (solid line) at 9 = offsets the positive contri- 
bution from Vt(r) (dashed line). The greater amplitude 
of Vt(r) in the bisymmetric model of NGC 2976 is there- 
fore due to the large non-circular motions in the inner 
parts that happen to be negative near the kinematic ma- 
jor axis because the m = 2 distortion is oriented close to 
this axis. 

Notice also that the V2,t and V<j,r components in Fig. [6] 
show, as they must ($2]), both ml = 1 and ml — 3 pe- 
riodicities, and that both are of similar amplitude (see 
also Fig. HJ. Yet their relative phases ensure that the 
net effect of the ml = 3 terms on Vmodei cancels almost 
exactly. A larger ml — 3 signal could arise if the V^t 
and V<2,r terms have different amplitudes, but their rel- 
ative phases always ensure at least partial cancellation 
regardless of the orientation of the projected bar. Thus 
one should not conclude that a very weak ml — 3 signal 
in the velocity map implies no significant bisymmetric 
distortion. 

5.2. Centrifugal balance? 

It is clear from Tableland Fig.[2]that both the bisym- 
metric and radial models are adequate parameterizations 
of the observed geometry and kinematics of NGC 2976. 
But does either model provide insight into its physical 
structure? 

The mean orbital speed, Vt (r) , in the radial model can 
balance the central attraction of the system only if either 
the non-circular motions are small, or V r (r) actually im- 



plies a real radial flow that somehow does not affect or- 
bital balance. The first possibility is not true, as we have 
confirmed (Figs. [3] and the large no n-circu lar motions 
found for r < 500 pc in NGC 2976 by ISBLBl . If V r (r) is 
attributed to radial flows that do not affect orbital bal- 
ance, then all of the detected gas in this quiescent system 
woul d be displaced on k pc scales in 1 — 3 Gyr; we agree 
with ISimon et al.l (2005) that this explanation is not vi- 
able. We thus conclude that although the optimal radial 
model is a reasonable statistical fit to the data and pro- 
vides strong evidence for non-circular motions, the fitted 
Vt(r) cannot be used to determine the mass distribution 
within NGC 2976. 

If the non-circular motions in NGC 2976 are domi- 
nated by an m = 2 perturbation to its potential, then 
the velocity profiles of the optimal bisymmetric model 
should better reflect the galaxy's structure than those 
of the radial model. While the fitted Vt(r) rises more 
steeply in the bisymmetric model, it is merely the aver- 
age azimuthal speed around a circle, not a precise indi- 
cator of circular orbital balance. It should be stressed 
that circles in the disk plane approximate streamlines 
only when non-circular motions are small. In a bar-like 
potential, the gas on the bar major axis will be moving 
more slowly than average, since it is about to plunge in 
towards the center, whereas gas at the same galactocen- 
tric radius on the bar minor axis will be moving faster 
than average, since it has arrived there from a larger ra- 
dius. Under these circumstances, it_ is not possible to 
assert that the azimuthal average, Vt, is exactly equal 
to the circular orbit speed in an equivalent azimuthally 
averaged mass distribution. The only reliable way to ex- 
tract the azimuthally averaged central attraction in this 
case is to find the non-axisymmetric model that yields 
a fluid dynamical flow pattern to match that observed, 
and to average afterwards. 

Despite these cautionary statements, we suspect that 
the Vt curve from the bisymmetric model provides a bet- 
ter estimate, than does that from the radial model, of the 
azimuthally averaged central attraction in NGC 2976. 

5.3. Evidence for a bar 

As discussed in §$T]&[3l the elliptical streams of fixed 
direction and phase in the bisymmetric model could be 
driven by either a triaxial halo or by a bar in the mass 
distribution. In either case, the distortion is significant 
only at r < 80" (1.4kpc; Fig.©) in NGC 2976, beyond 
which the flow appears to be near circular. 

The aspherical halo interpretation therefore requires 
the halo that hosts NGC 2976 to have an aspheric- 
ity that inc reases for decr e asing r. Such an idea was 
roposed bvlHavashi et all (120071), al t hough other work 



proposed by riavasln et al.l \\Z\J 
(iDubinskil 119941: iGnedin et al. 



iBerentzen fc Shlosmanl 2006 : iGustafsson et al 



2004 : iKaz antzidi s et al 



2006J) has indicated a tendency for disk assembly to cir- 
cularize the potential. 

We therefore favor the interpretation th a t NGC 2976 
hosts a bar. Menen dez-Delmestre et all ((2007) have 
examined the Two Micron All Sky Survey (2MASS: 
iSkrutskie et "all[200l J, H and K s images of NGC 2976 
to search for a bar. Their fits to this photometry reveal a 
radi al variation in el lipticity of amplitude Ae > 0.1 (see 
also lSimon et al.ll2005f ). and their visual inspection of the 
images reveals a "candidate" bar with PAb ar = —43° and 
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TABLE 1 
Minimization Results 



Model e d <j>' d x c y c V sys <j>' b X ? min " (|AV n |) 

(°) (") (") (fans- 1 ) (°) (kms" 1 ) 

(1) (2) (3) (4) (5) ' (6) (7) (8) (9) (10) 



radial 0.568 ±0.007 -36.0 ± 0.6 -1.7 ±0.3 0.8 ± 0.4 0.3 ± 0.4 ■■• 1.35 1034 3.4 

bisymmetric 0.556 ±0.007 -37.6 ± 0.6 -1.9 ±0.3 1.2 ± 0.3 0.7 ±0.4 -45 ± 4 1.20 1009 3.1 



Note. — Col. (1): Model. Col. (2): Disk ellipticity. Col. (3): Disk PA, measured North East to the receding 
side of the disk. Col. (4): Right ascension of disk center, relative to the phot ometric center 09 h 47 m 15? 3 iSBLBT) , Col. 
(5): Declination of disk center, relative to photometric center 67° 55' 00'.' 4 (SBLB). Col. (6): Disk systemic velocity, 
heliocentric optical definition. Col. (7): Bisymmetric distortion PA, measured North — > East. Col. (8): Minimum 
value of x? (eq-EJ obtained. Col. (9): Number of degrees of freedom in the minimization. Col. (10): Amplitude of the 
average (data - model) residual. 



semi-major axis a = 72 ± 5" (see their table 2). Their 
estimated PAbar is fully consistent with our kinematic es- 
timate <j)' b (Table [T|) , and a compares well with the range 
of r where V2,*(r) and V2, r (r) are non-zero (Fig. [3^). Fur- 
thermore, PAbar arid <f>' b are roughly coincident with the 
apparent ma jor axis of the CO distribution in NGC 2976 
(see fig. 4 of ISBLBl - ). which suggests that the molecular 
gas density is larger along this PA than elsewhere in the 
disk. Thus the 2MASS photometry and CO morphol- 
ogy provide strong supporting evidence that NGC 2976 
contains a bar with the properties implied by our bisym- 
metric model. 

5.4. Mass components 

Our fits have revealed strong non-circular motions in 
NGC 2976 that appear to result from forcing by a bar. 
While Vt (r) in the bisymmetric model better reflects the 
azimuthally averaged mass distribution than its counter- 
part in the radial model, precise statements about the 
mass budget in NGC 2976 are hampered by our lack of a 
reliable estimate of the circular orbital speed curve (see 

The amplitude of the non-circular motions in 
NGC 2976 implies a relatively large bar mass, which 
in turn suggests that the disk itself contributes signifi- 
cantly to the central attaction. It is therefore likely that 
the baryons in NGC 2976 dominate its kine matics well 
beyond the r ~ 500 pc suggested by the fits of lSBLBl . In- 
deed, the steeper rise of Vt( r) in th e bisymmetric model 
relative to that deduced by ISBLBl would allow a larger 
disk mass-to-light ratio (Ai/L) to be tolerated by the 
kinematics. This conclusion eases the tension between 
their dynamical upper bound on the stellar M/L and 
that expected from stellar population s ynthesi s for the 
observed broadband colors (see §3.1.1 of SBLBl). 

We defer the detailed mass modeling of NGC 2976 re- 
quired for quantitative estimates of its mass budget to a 
future paper. Such a study would be assisted by addi- 
tional kinematic data from extant H I aperture synthesis 
observations, as well as by decompositions of publicly 
available infrared photo metry from the Spitzer Infrared 
Nearby Galaxies Survey (|Kennicutt et al1l2003h . 

5.5. Other galaxies 

We suggest that our approach could be useful for char- 
acterizing the non-circular motions detected in other 
galaxies, particularly in low-mass syste ms where the 
report ed non-circular m otions are large (|Swaters et al.l 
20031 iSimon et al.l[2005t iGentile et al.ll2007h . It is more 



direct than interpretations of velocity field Fourier coef- 
ficients in the weak perturbation limit, yields physically 
meaningful kinematic components for systems with bar- 
like or oval distortions to the potential, and its applica- 
tion is much simpler than that of a full fluid-dynamical 
model. 

We have shown that the velocity field of NGC 2976, 
when fitted by_our bisymmetric model, reveals a steeper 
inner rise in Vt(r) than in previous analyses by other 
methods. Similar findings have been reported by Hayashi 
& Navarro (2006) and Valenzuela et al. (2007) for other 
systems. In NGC 2976, the reason for this difference 
(Fig. [6]) is that the V^t terms happen to partly cancel 
the V t terms, because they have opposite signs on the 
projected major axis when the bar is oriented close to 
this direction. It should be clear, however, that the Vi,t 
terms will have the opposite effect if the bar is more 
nearly aligned with the projected minor axis. Thus even 
if the non-circular motions detected in other systems re- 
sult from bars in the potential, it is unlikely that our 
bisymmetric model will always cause Vt(r) to rise more 
steeply than found previously. In any event, it should be 
clear that when large non-circular flows are present, the 
mean orbital speed derived from models that use epicy- 
cle theory can yield a very misleading estimate of the 
interior mass needed for centrifugal balance. 

6. CONCLUSIONS 

We have presented a new method for fitting 2-D ve- 
locity maps of spiral galaxies that are characterized by 
non-circular motions. We suppose the potential to con- 
tain a bar-like or oval distortion that drives the gas in 
the disk plane on an elliptical flow pattern of fixed ori- 
entation, such as could arise from a triaxial halo or a bar 
in the mass distribution of the baryons. Our model has 
important advantages over previous approaches since it 
is not restricted to small non-circular motions, as is re- 
quired when epicycle theory is employed, and we do note 
invoke large radial flows that have no clear physical origin 
or interpretation. 

Our bisymmetric flow model can be fitted to data 
by a generalizatio n of the technique developed by 
iBarnes fc Sellwoodl (|2003h . The fit extracts multiple 
non-parametric velocity profiles from an observed veloc- 
ity field, and we employ a bootstrap method to estimate 
uncertainties. 

As an example, we have applied our technique to the 
Ha and CO kinematics of NGC 2976 presented bv lSBLBl . 
We show that the bisymmetric model fits the data at 
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Fig. 6. — Projected contributions from different kinematic components at r = 20" (345 pc) in the optimal (a) radial and (b) bisymmetric 
models, fn [BJi, the dashed line shows the angular dependence of the projected rotational velocity term in the radial model relative to the 
kinematic major axis (2 nd on the right-hand side (RHS) of eq. [7J, and the dash-dotted line shows that of the radial velocity term (3 rd on 
the RHS of eq.0. In[Bf), the angular dependence of the components in the bisymmetric model are plotted relative to <j>' d along the bottom 
horizontal axis and <f>' b along the top horizontal axis. The dashed line shows the rotational velocity term (2 nd on the RHS of eq. |5j, and the 
dash-dotted and dash-dot-dotted lines show the tangential and radial bisymmetric terms, respectively (3 rd and 4 th on the RHS of eq. [5}. 
The solid lines in both panels shows the net projected model velocity relative to V sys . 



least as well as the ringfit procedure implemented by 
these authors that invokes large radial velocities, which 
we are also able to reproduce by our methods. Both 
the bisymmetric and radial models reveal large non- 
circular motions in NGC 2976, but the derived mean 
orbital speed profiles V t (r) differ markedly between the 
two cases. We explain the reason for this large difference 
in H 

When disks are observed in projection, kinematic dis- 
tortions with intrinsic sectoral harmonic m cause az- 
imuthal variations of orders ml = m ± 1 in line-of-sight 
velocity maps. Our analysis of NGC 2976 clearly demon- 
strates that m! = 1 distortions to its velocity field can 
be fitted by a bisymmetric distortion to the potential, 
which we regard as more physically reasonable than ra- 
dial flows. We show in Fig. [6] that m' = 3 distortions 
should be small in the bisymmetric model; this is be- 
cause the m! = 3 variations in the radial and tangen- 
tial components project out of phase. They will cancel 
exactly only when of equal amplitude, which should be 
approximately true in the rising part of Vt(r). 

We suggest that NGC 2976 hosts a strong bar oriented 
at ~ 17° to the projected major axis. Our int erpretation 
is supported by its CO morp holoev (ISBLBT) and more 
strong ly by the results of iMenendez-Delmestre et all 
(|2007h . who analyzed the 2MASS photometry of 
NGC 2976 and found a bar whose size and orientation 
are similar to those required by our bisymmetric model. 

We find that the mean orbital speed in NGC 2976 
rises more steeply than i ndicated by previous studies 
(jSBLBt ISimon et al.ll2005h . While V t (r) in our bisym- 
metric model is not an exact measure of the circular or- 
bital speed in the equivalent axially symmetrized galaxy, 
we regard it as a better approximation to this quan- 
tity. Since the strongly non-circular flow pattern implies 
a massive bar, which in turn suggests a massive disk, 
we exp ect a larger baryonic mass than was estimated by 
SBLB. It is likely, therefore, that most of the increased 
central attraction required by our more steeply rising 
Vt (r) will not reflect a corresponding increase in the den- 
sity of the inner dark matter halo, but will rather ease 
the tension between maximum disk fits to its kinemat- 



ics and M/L predictions from broadband photometry 
(SBLB). Indeed, since non-circular motions are detected 
throughout the region r < 80" (1.4 kpc), it seems likely 
that the luminous matter in NGC 2976 is an important 
contributor to the central attraction at least as far out 
as this radius. Detailed mass models of this system are 
forthcoming. 

Application of our method to other galaxies will not 
always result in a steeper inner rise in the mean orbital 
speed. We find this behavior in NGC 2976 only because 
the bar is oriented near to the projected major axis. Ne- 
glect of non-cirular motions, or application of a radial 
flow model, when the bar is oriented close to the pro- 
jected minor axis will lead to an erroneously steep rise in 
the apparent inferred mean orbital speed, which will rise 
less steeply when our model is applied. 

We thank Josh Simon for providing the data for 
NGC 2976, and Alberto Bolatto for help in interpreting 
the measurement uncertainties. We also thank Alberto 
Bolatto and Josh Simon for helpful comments on the 
manuscript. KS is a Jansky Fellow of the National Ra- 
dio Astronomy Observatory. JAS is partially supported 
by grants from the NSF (AST-0507323) and from NASA 
(NNG05GC29G). 
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APPENDIX 

KINEMATIC MODEL WEIGHTS 

We use the same notation and parameter definitions as in fJ3]and Fig. [TJ Let (x n ,y n ) be the location, in the sky 
plane, of the nth measured velocity D n , where x points West and y points North. Let <f>' d be the PA of the projected 
major axis. The coordinates (x e ,y e ) centered on {x c ,y c ) and aligned with the projection are: 

x e = -(x n - x c ) sin0' d + (y n - y c ) cos cj)' d 

y e = -(x n - x c ) coscj)' d - (y n - y c ) sin^. (Al) 
The radius r n of the circle in the disk plane that passes through the projected position of (x n ,y n ) is given by 

+ • (A2) 

where is the ellipticity of the circle caused by projection; i.e. — 1 — cos i for the thin gas layer considered here, 
where i is the galaxy inclination. Thus the PA, 0, of the measurement at (x n ,y n ) relative to the disk major axis (see 
Fig. [J) satisfies 

cos0= — , (A3) 
r n 

Ve 



(1 - ed)r v 



(A4) 



We tabulate each of the non-parametric velocity profiles at a set of radii in the disk plane that project to ellipses 
on the sky with semi-major axes {ak'm}- The elements of {14} include all K — J2m K'm values from the M velocity 
profiles combined. In principle, we could employ different numbers of rings K' M , with differing choices for the {ak'm} 
for each velocity profile, but here we evaluate all profiles at the same K' locations {aw}- There are K [ radii that 
define Vt(r) in our models (cqs. [5]and[7|), but we restrict the number of rings describing the non-circular components 
to K" < K' when the latter are not well-constrained in the outer disk. In addition, we include the systemic velocity 
V sys a s the K th element of {14}, with u>K, n = 1 f° r all D n . 

We use linear interpolation between the two rings that straddle each data point. If < r n < dk+i, the values of 
{wk, n } for that velocity profile and that D n are 

W k ,n = U 7 W n 

V oa k J 

( r n ~a k \ 
Wk+X,n = [—fa J W n 

Wfc',n = for k'^k, k+1 , (A5) 



where Sak = flfc+i — afc is the ring spacing. The parameter W n in eq. IA5I is the combination of trigonometric factors 
in the dependence of {wk, n } on the projection geometry. It is clear from eqs. [5] and [7] that we require a different W n 
for each D n and for each velocity profile in the bisymmetric and radial models. 
For the Vt component in both models, we have 

W n = sin i cos 6 , (A6) 
where cosO is given in eg. IA3I for the V% t component in the bisymmetric model 

W n = sin i cos(26» b ) cos , (A7) 

and for the Vi r component 

W n = sin i sin(26» & ) sin 6 . (A8) 
Finally, W n for the V r component in the radial model is 

W n = sin i sin . (A9) 
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